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Abstract 



The exact Green's functions of the periodic Anderson model for U ^ oo 
are formally expressed within the cumulant expansion in terms of an effective 
cumulant. Here we resort to a calculation in which this quantity is approxi- 
mated by the value it takes for the exactly soluble atomic limit of the same 
model. In the Kondo region a spectral density is obtained that shows near 
the Fermi surface a structure with the properties of the Kondo peak. Approx- 
imate expressions are obtained for the static conductivity cr(T) and magnetic 
susceptibility x(^) of the PAM, and they are employed to fit the experimental 
values of FeSi, a compound that behaves like a Kondo insulator with both 
quantities vanishing rapidly for T ^ 0. Assuming that the system is in the 
intermediate valence region, it was possible to find good agreement between 
theory and experiment for these two properties by employing the same set 
of parameters. It is shown that in the present model the hybridization is 
responsible for the relaxation mechanism of the conduction electrons. 
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I. INTRODUCTION 



In the present work we discuss approximate Green's Functions (GF) of the Periodic 
Anderson Model (PAM) that use the atomic hmit as a starting point. We employ these GF 
to calculate the static magnetic susceptibility and the resistivity of FeSi. 

FeSi has rather unusual magnetic properties,§ and in particular a static susceptibility 
x{T) that has a maximum x{Tm) at about = 536 K and vanishes for T — when 
the low temperature Curie tail is subtracted.B Several models were studied by Jacarino 
et.al.,0 and they found that a simple one, with two very narrow rectangular bands separated 
by a gap, could be used to fit x(T). The properties of FeSi are very similar to those of 
the Kondo insulators,iHi and a very simple model that describes most of their properties 
would consist of two hybridized bands with two electrons,! i.e.: an intrinsic semiconductor 
with an hybridization gap. Both x{T) and the resistivity p{T) of FeSi have activation laws 
with characteristic energies of order 0.1 eV, and although band structure calculationJl give 
comparable semiconducting gaps they can not reproduce the large values of x{T)- Also the 
measurements of infrared and optical reflectivity can not be described with the predictions 
of the band calculations,! and the PAM seems a more adequate model, both because it 
becomes the previous model when the Coulomb repulsion between the localized electrons is 
neglected, i.e.: U = 0, and because the correlations present in the PAM could explain some 
of the observed features of FeSi. 

A two-band Hubbard model! has been already employed to study FeSi, and it is equal to 
the PAM when the dispersion in one band and the Coulomb repulsion in the other are zeroed. 
Another variant of the PAM was to add dispersion to the band of localized electrons, and 
this model has been studied both with U — > oM and with finite U.i It has been suggested! 
that it would be interesting to describe the Kondo insulators employing the PAM with 
U oo, and this approach is presented here to calculate both x{T) and p(T) employing 
approximate! GF. This model has been also used to study the Kondo insulator CeaBi4Pt3, 
employing the slave boson technique in the mean field approximation.li3liil Varmali^ argues 
that the strict Kondo lattice (i.e. with integer f occupation and therefore without charge 
fluctuations) is inappropriate for these systems, because it would be rather unlikely to find 
the chemical potential just in the hybridization gap, and that one would more likely find 
this situation in a mixed valence system. Our calculation shows that with the simplified 
model presented here, it is possible to give a fair description of the x{T) and p(T) of FeSi 
in a typical intermediate valence situation.El 

In the study of solid state systems it is sometimes interesting to focus on the local states 
of the ions placed at the different sites of the crystafl The space of the local or ionic states 
associated to a given site contains many states that are usually of little interest, generally 
because their occupation can be neglected at fairly low temperatures or frequencies when 
they are too far apart in energy from the ground state. It is then useful to eliminate these 
states from the model of the system, and the Hubbard operatorsBlll are very convenient 
for that purpose. In the Anderson lattice we have a broad band of conduction electrons, 
identified by a subindex c, and the local d states, which will be identified with the subindex 
/ for convenience. At each site j of the lattice there are four local states: the vacuum state 
I j, 0), the two states | j, a) of one electron with spin component a and the state | j, 2) with 
two local electrons. When U ^ co the state | j, 2) is empty, and one can use the Hubbard 
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operators to project it out from the space of local states at site j. One difficulty is that 
the usual expansions employed with the Fermi or Bose operators are not valid for these 
operators, and Hubbard introduced for his model of correlated electron^ a diagrammatic 
expansion with cumulant^l which uses the electron hopping as a perturbation, and becomes 
the usual expansion when U = 0. The GFs employed in the present work are based on an 
extensionS of Hubbard's cumulant expansion that is valid for the Anderson lattice and uses 
the hybridization as perturbation; employing this expansion it is possible to express the 
exact GF in terms of an unknown effective cumulant&El M2^J{uj). In this work we shall 
approximate the effective cumulant by M^^^lz), which is obtained from the exact solution of 
the Anderson lattice in the atomic limit, namely when the band of uncorrelated electrons has 
zero width!. The spectral density obtained in this approximation, for typical values of the 
system parameters, shows a structure close to the chemical potential /i that corresponds to 
the Kondo resonance and affects the physical properties at low temperatures. This structure 
was absent from the GF derived from the cumulant expansion when only cumulants up to 
fourth order were employedlllS, and this was the motivation to approximate the effective 
cumulant by a method that would include all the higher order cumulants that were absent 
in our previous calculations. Because of its atomic character, the approximate effective 
cumulant in our method is independent of the wave vector. Our GFs are therefore closely 
related to those employed by the dynamical mean field theoryl^ of the infinite dimensional 
problem, but rather than using a self consistency condition we use physical considerations 
to chose the parameters that define the effective cumulant. The dynamical mean field theory 
has been also employed to study the magnetism of the PAMcil as well as the transfer of 
spectral weight in the spectroscopy of correlated electron systems described by this model.E^ 

Both the dynamical susceptibility and the conductivity require two-particle GFs in the 
theory of linear response, but we shall obtain approximate expressions of the static quantities 
employing the one-particle GFs introduced in this work. With these GF we obtain the total 
number of spin up and of spin down electrons in the presence of a weak magnetic field, and the 
static susceptibility x{T) is then proportional to their difference divided into the magnetic 
field; we can then compare the ratio x{T)/x{Tm) with the corresponding experimental value. 
To simplify the calculation of the effective cumulant we assume equal gyro magnetic factors Qf 
and Qc, although the extension to different gj and Qc would not present essential difficulties. 

It has been shownil that the limit d = oo can be used to give an approximate description 
of three-dimensional systems, and we shall then use an expression of a (T) that is valid in 
infinite dimension in our calculation for FeSi. This expression is derived from the well known 
Kubo formula, the vertex corrections cancel out when d = oo,0 and only the one-particle 
Gc^aO^, z) are then necessary to calculate a (T). In this expression there are explicit sums over 
k, but when nearest-neighbor hopping in a simple cubic lattice is considered, it is possible to 
derive expressionJUil that depend on k only through the unperturbed conduction electron 
energies e (k). As a further simplification that would not change the results in an essential 
way, we shall use a rectangular band with —W < e (k) < W . 

Using the expressions of a (T) and xiT) discussed above we fitted the experimental 
magnetic susceptibilityS x{T)/x{Tm) and the static resistivityi p (T) = 1/cr (T) of FeSi with 
the same set of parameters in a typical situation of intermediate valence, obtaining a fairly 
reasonable agreement with the experimental values. To adjust the x{T) at high T, it was 
necessary to assume that the thermal expansion affects the value of the system's parameters. 
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In Section |I| we discuss the Periodic Anderson Model (PAM) and the approximate one- 
particle GF employed in this work. In Section III we analyze the static magnetic suscep- 
tibility x{T) and the static resistivity p (T) of FeSi. Conclusions are presented in Section 

m 



II. GREEN'S FUNCTIONS FOR THE PERIODIC ANDERSON MODEL 

As discussed in Section | we employ Hubbard's operators.0 In the general case there is 
a fixed number n of orthogonal states { \j, a) } (identified by indices a) that span a space 
Sj^n at each site j, where j = 1,2, . . . , Ng and A^^ is the number of sites. To each site j we 
associate the Hubbard operators 

Xj,ab = \j,a){j,b\ , (2.1) 

which transform the state \j,b) into the state \j,a) , i.e. Xj^ab\j,b) = \j,a)- The product 
rules for two operators at the same site are given by 

Xj,ab Xj cd = Sb,c Xj^ad , (2.2) 

and we chose properties equivalent to those of the usual Fermi or Bose operators when they 
are at different sites. We then say that Xj^at is of the "Fermi type" ( "Bose type" ) when the 
number of electrons in the two states | j, a) and | j, b) differ by an odd (even) number. For 
j 7^ j' we then use {Xj^ab, Xji^cd} = when the two operators are of the "Fermi type" and 
[Xj^ab, Xj/^cd] = when at least one is of the "Bose type" (as usuaS [a,b] = ab — ba and 
{a, b} = ab + ba). 



A. The Anderson lattice for [/ — > oo 

The Anderson lattice with finite U is described by the Hamiltonian 

-f^ = X! -^k,CT C'LCkCT + fjafja + X! ^ fjafj^rfjafja + 

kcr j(7 j 

E + H.C.) , (2.3) 

jka 

where Cl^ {Cka) is the usual creation (annihilation) operator of conduction band electrons 
with wavevector k and spin component ah/2, where a = ±1. The /j^ fj<T correspond 
to the local f or d electrons at site j , and 

^,,k,. = V^(A:)exp(^k.Rj) . (2.4) 

When f/ = Eg. ( |2.3| ) describes two hybridized bands of uncorrelated electrons. 

The state space of the f-electrons at each site j is spanned by the four states] j, 0), | j, a) 
and I j, 2), with a = ±1. The state | j, 2) is empty when U oo, and we shall consider a 
reduced space of states by projecting | j, 2) out, so that £j^n is a three dimensional space. 
To make the connection with the Hubbard operators one could substitute the identityB^ 
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fja — ^j,Oa + O'-^j.CTci (2-5) 

into Eq. (^.31), where the factor a is necessary to obtain the correct phase of the states. 
Ehminating the Xj^-^d and Xj^22 one obtains the projection of H into the reduced space, 
namely 

kcr J.CT 

XI (^i,k,a -^],0a C'k.a + Vj*^^^ C'k.a ^i,Oa) ■ (2.6) 
jkcr 



1. The cumulant expansion 

The cumulant exnansion has been employed by several authors to study Ising's and 
Heisenberg's models,l23 while Hubbardli extended the method to a quantum problem with 
fermions. In this technique the cumulant averages^ are used to rearrange the usual perturba- 
tive expansion,0 and it is possible to derive a diagrammatic expansion involving unrestricted 
lattice sums of connected diagrams, that satisfies a linked cluster theorem. This technique 
was extended to the Anderson lattice,^ and a brief description is given here. The method 
employs the Grand Canonical Ensemble of electrons, and it is then convenient to introduce 

where is the chemical potential and Va is the number of electrons in the state | j, a) . The last 
term in Eq. (|2.6| ) will be considered as the perturbation, and the exact and unperturbed 
averages of any operator A are respectively denoted hy < A >t-i and < A >. It is also 
convenient to introduce 

Ej^a = Ej^a - fJ'^a (2.8) 

ei,a = E]^^- fi , (2.9) 

because these are the forms that consistently appear in the calculations. 

The Matsubara expansion^ is employed, so that r is an imaginary time in the GFs 

(Xj air) Xji a'('r')) ) , where a = (a, b) identifies the transition b a, and 

Xj^air) = exp (tH) X^- „ exp {-tH) . (2.10) 

The subindex + in the definition of the GF indicates that the operators inside the parenthesis 
are taken in the order of increasing r to the left, with a change of sign when two Fermi-type 
operators have to be exchanged to obtain this ordering. The inverse of Plank's h and of 
Boltzmann's /c^ are usually included into the real and imaginary temperatures T and r, 
as well as in several other parameters, so that all of them are given in terms of a common 
energy unit. 
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Some of the infinite diagrams which contribute to the GF ^(^Xj Q,(r) Xj/ a/(r')j^y are 

shown in figure |l|, and the full circles (f-vertices) correspond to the cumulants of the f- 
electrons. Each line reaching a vertex is associated to one of the X operators of the cumulant, 
and the free lines (i.e. those that do not join an empty circle) correspond to the external X 
operators appearing in the exact GF. An explicit definition of the cumulants can be found 



in the references p!7| , p!8| , |34| , and they can be calculated by employing a generalized Wick's 
theorem.tJES 

The first diagram in figure |l|a corresponds to the simplest free propagator 
(Xj Q,(r) Xj/ a')^^, and the second diagram in that figure has an empty circle (c- 
vertex) that corresponds to the conduction electron cumulant, equal to the free 

propagators^ (^Cfco-(T) C^^^^^ = G'°^(k,r). The interaction is represented by the lines (edges) 

joining two vertices and, because of the structure of the hybridization, they always join a 
c-vertex to an f- vertex; the number of edges in a diagram gives its order in the perturbation 
expansion. 

Cumulants containing statistically independent operators are zero, and those appearing 
in the present formalism (with the hybridization as perturbation) vanish unless they contain 
only X operators at the same site or only C or C"'" operators with the same k and a. The 
only non-zero c-cumulants are of second order, because the uncorrelated c-operators satisfy 
Wick's theorem. On the other hand, the f-vertices can have many legs, all corresponding 
to X operators at the same site, like the fourth and sixth order cumulants appearing in the 
rather more complicated diagram shown in figure |I]c. 

All the infinite diagrams that contribute to the GF with cumulants of at most second 
order are shown in figure |T|a, and this family is the "chain approximation" (CHA), which 
gives the exact solution of Eq. ( p.3| ) when there is no Coulomb correlation {U = 0). When 
the spin is eliminated from the problem, the Hamiltonian of equation ( p. 61) corresponds to 
a system of two hybridized bands without any Coulomb repulsion (there can be only one 
or zero f-electrons at each site), and the CHA is again an exact solution.lli The diagrams of 
the c-electron propagator in the same approximation are shown in figure 

In the Feynmann perturbation expansion. Wick's theorem is valid and only second order 
propagators appear, while the interactions are provided by the Coulomb interaction. In the 
present treatment, the U disappears in the limit f/ ^ oo (or is included in the unperturbed 
Hamiltonian when U is finite), and the correlations appear through the non-zero cumulants 
of X operators with order greater than two, which include propagators of two or more 
particles. In the Feynmann expansion of the one particle GF, the two particle GF appear 
in the self-energy, which contains all the correlations. 

2. The spectral density of the GF and the occupation number 
In the Anderson lattice with U ^ oo one can introduce one-particle GFs of local electrons 

((X,,o.(r) X,,,o(r')) + )^ , (2.11) 
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as well as GFs for the c-electrons {^{Ckair) CL(r'))^^^ and "crossed" GFs of the type 
(^Xj Oo-(T) C|,^(r')^ \ , all of them defined in the intervals < r, r' < /? = 1/T. It is 



n 

possible to associate a Fourier series to these GFs because of their boundary condition in 
this variable,lli and the coefficients correspond to the Matsubara frequencies Wj, = Tiiy/P 
(where u are all the positive and negative odd integer numbers). One can also transform 
the GF to reciprocal spaced so that 

(X„o.(r) X,, .o(r'))+)^ = ^ E E ^xp [z (k.R, - k'.R; 



T + t')] (^(Xk,o<7(u;i.) -^k',<7o(t^i.'))+/^ 5 (2-12) 
and because of the invariance against time and lattice translation^lil 

{{X\ifla{^u) -^k',(7o(t^i^')) + )^ = ^//,o-(k, ^^i/) 5k',k 5v+v',0 , (2-13) 

where 5k',k and 5^+i.',o are Kronecker's deltas, because the k and Ui, are discrete variables. 
Transforming the eigenstates of the c-electrons to the Wannier representation, one obtains 
the equivalent relations for Gcc,cr(k, u^) and G/c,o-(k, uji^). Considering that the coefficients of 
the r Fourier series for each k are the values of a function of the complex variable z = u + iy 
at the points = i u^, it is possible to make the analytic continuation to the upper and 
lower half-planes of z in the usual way,@ obtaining, e.g. from the G ff^cJ^,i^u), a function 
Gff^aO^, z) which is minus the Fourier transform of the double time GF.E3 

If we assume the system to be uniform, the occupation of the local state {Xj^^^) does 
not depend on j, and it is given by 



^/,'^ = / Pffia{i^) fri^) duj , (2.14) 

where 

fT{z) = [l + expiP z)]-' (2.15) 

is the Fermi function and 



- hm Im\^J2GffA^,oo + t\r^\)] (2.16) 



is the spectral density associated to the transition a ^ 0, abbreviated with (0, a). Using 
the same pf^o„{u) it is also possible to obtain the occupation (Xj^q) = Uf^ of the empty 
state I 0), namely 

/oo 
pffiaii^) (1 - fAi^)) duj . (2.17) 
-oo 

All the corresponding quantities for the conduction electrons are obtained by replacing 
Gcc,a{K z) for G'//,CT(k, z). 
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The f-electron GF is given in the CHA (cf. the diagrams in figure |l|a) by 



where the energies £:i,o-(k) and £2,0- (k) are the two elementary excitations with wave vector 
k and spin component cr, resulting from the hybridization of a band £k(T and a dispersionless 
band of energy ej = Ej„ — ft, with a reduced hybridization constant ^/D^V{k). These 
energies are given by the two roots of {z — Sf){z — 6]^^) — D„ \V{k)\^ = 0, where 

/oo 
Pf,Oa{u;) du = {Xoo + X^„) . (2.19) 
-00 

The spectral density for the unperturbed f-electrons is a 5-function at Sf, and pf^Q^ioj) 
becomes two bands with a gap centered at and roughly proportional to the reduced 
hybridization constant \/D„V in the CHA. When the system is in the "Kondo region" in 
which the local state has nearly the maximum occupation compatible with the paramagnetic 
state, viz. n/^o- = 0.5 in the average, a narrow temperature-dependent peak, the "Kondo 
peak", should appear very close to the chemical potential, i.e. near to = in the variables 
we use. This peak is responsible for many properties of the "heavy fermions",0 and it is 
closely related to the increase in resistivity when T decreases, observed in many metals 
with magnetic impurities at low T. In combination with the effect of phonons, that make 
the resistivity increase with T, it is responsible for the minimum in resistivity, named the 
Kondo effect, that is observed in those systems. To explain this behavior it was essential to 
consider a third order perturbation that includes spin-flip processes.0 These processes are 
absent from the CHA, because diagrams with only second order cumulants do not include 
them, so that the absence of the Kondo peak in that approximation is not surprising. The 
GFs proposed in this work avoid the laborious explicit calculation of higher order cumulants 
by including all of them in an approximate way. 



3. The exact Green's functions 

In the calculation with the usual Fermi or Bose operators, the one-particle propagator 
of the f-electron is given by a sum of diagrams of the same type0 shown in figure [l|a but 
with each vertex corresponding to the sum of all "proper" (or irreducible) diagrams.0ii 
The same result is found in the cumulant expansion of the Hubbard model for d cx),Ei@ 
when the electron hopping is employed as perturbation. The vertices then represent an 
"effective cumulant" M2^J [z), that is independent of k because only diagrams of a special 
type contribute to this quantity for d — > cxd. 

In the cumulant expansion of the Anderson latticelll we employ the hybridization rather 
than the hopping as a perturbation, and the exact solution of the conduction electrons 
problem in the absence of hybridization is part of the zeroth order Hamiltonian. For this 
reason it became necessary to extend Metzner's derivational to the Anderson lattice, and we 
have shownc3 that the same type of results obtained by Metzner are also valid for this model. 
These results had been used! to obtain the exact GF employed in the present work, but since 
then we realized that the expression of the exact GF is valid for all dimensions and it is not 
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necessary to assume infinite dimension in that part of the derivation. As with the Feynmann 
diagrams, one can rearrange all those that contribute to the exact Gff^„(k, z) by defining an 
effective cumulant M2^/ (k, z), that is given by all the diagrams of Gff^(j(k, z) that can not 
be separated by cutting a single edge (usually called "proper" or "irreducible" diagrams). 
The exact GF Gjj o-(k, z) is then given by the family of diagrams in figure |l|a, but with the 
effective cumulant M2^/ (}i,z) in place of the bare cumulant M2^{z) = —D'^/{z — ej) at all 
the filled vertices. The exact GF for the f electron is then written as 

GffA^, z) = M|,{^(k,.) 2 ro\v , M^fffv , ' (2.20) 

1- I V{k) |2 G'°^(k,z) M2-;J{k,z) 

where G°^^(k,z) = —l/{z — £(k)) is the frequency Fourier transform of G°^(k, r), and in a 
similar way one obtains the exact GF for the c-electron, namely 

G^^Jk,z) = — Tf . (2.21) 

z-eik)+\Vik)\^ M^fJ{k,z) ^ ' 

It is clear that for arbitrary dimension we have not gained much with Eq. ( p.20| ) , because 
the calculation of M2^/{k,z) is as difficult as that of G'// a-(k, z). We circumvent this obstacle 
by resorting to replace M2^/ {k,z) by the corresponding quantity M2^^{z) of an exactly soluble 
model, which is the atomic limit of the Anderson lattice. The hopping is neglected in this 
system, described by the Hamiltonian of Eqs. ( |2.3| ) or ( ^.6] ) with _Ek o. = Eq, and it already 
contains the basic physics of the formation of the singlet ground state and of the appearance 
of the Kondo peak in the PAM, as it is discussed in review articles.^ Because of its atomic 
character, the approximate effective cumulant M2^^{z) thus obtained is independent of k, 
and can be calculated exactly as shown in the next section. Being a special case of the PAM, 
it implicitly contains all the higher order cumulants. 

The effective cumulant is also independent of k when d ^ oo, and the very successful 
"dynamical mean field theory" also employs GFs corresponding to those in Eqs. ( [^.20| , py^TD . 

B. The exact solution in the atomic Umit 

The exact solution of the local problem has been already used in different ways to 
studv the Anderson lattice. The limit [/ — >■ oo was studied in the intermediate valence 
casetZl by considering only the lowest four eigenstates of the local Hamiltonian and the mag- 
netic instabilities and susceptibility were discussed employing the resulting self-consistent 
Hamiltonian.^ Considering only the atomic limit, viz. taking E]^^ = E^, Alascio et.alE^S 
studied the model for the whole range of parameters, showing that "most of the essen- 
tial characteristics" of these systems "are present in this crudely simplified Hamiltonian" . 
Simoes et.al.0 employed the atomic limit together with a diagrammatic method,ii that is 
essentially equivalent to our CHA, considering both the hopping and the hybridization as 
perturbations. An important improvement of the technique was to apply the same diagram- 
matic expansion to the exact solution of the atomic limit J^H employing only the hopping 
as Derturbation, and this technique has also been applied to study the problem with finite 
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As discussed in the previous section we introduce the exact expression for the GF that 
is given in terms of an effective cumulant M2'^J (k, z), however replacing this quantity by the 
approximate M2*^{z). As this treatment was derived from the diagrammatic expansionlll 
which uses the hybridization as perturbation and employs the exact solution of the uncor- 
related conduction band, it seems a better starting point than considering the hopping as a 
perturbation, because the hybridization is usually rather smaller than the bandwidth. The 
atomic limit has also been employed, within the framework of the dynamical mean field 
theory, to study the transport properties of the symmetric PAM.0 



1. The exact GF in the atomic limit 

Taking i^k.o- = Eq and introducing a local hybridization l^-,k,o- = the eigenvalue 
problem of Eqs. (p.3|) or ( p.6|) has an exact solution,^ and the GFs can be calculated 
analytically. As the problem is fully local, one can use the Wannier representation for the 
creation and annihilation operators Cj^ and o- of the c-electrons, and write = J2j Hj, 
where Hj is the local Hamiltonian 

Hj = ^ |-Eo C'j>Cj,<T + Ej^^ Xj^^a + (Vj^a xIq^ Cj^a + V*^ Cj^^ ^jfia)] , (2.22) 

cr 

and the sub index j can be dropped because we assume a uniform system. 

We shall denote with | n, r) the eigenstates of the Hamiltonian Hj with eigenvalues En^r, 
where n is the total number of electrons in that state, and r characterizes the different states. 
These eigenstates satisfy 

H \n,r) = En^r \ 'n^r) , (2.23) 

where H. is given in Eq. ( |2.7| ) and Sn^r = -^n.r — ""-/^ (cf- Eq. ( |2.8| )). In Table |I| we give the 
properties of the | n^r) states: number of electrons ra, name of the state r, z component of 
spin and en,r = En,r — nfi. The twelve eigenvalues en,r of Hj are represented in figure 0, 
and those corresponding to different occupations n = 0, 1, 2, 3 are drawn in different columns. 
The states are identified in the figure by the numbers r above the levels, and the lines joining 
different levels correspond to the possible transitions that contribute to the GF. 

It is now straightforward to express the Fourier transform of the f-electron GF in the 
form 

^(X,,o.(o;,) XlM))\ = ^{u^s+uj's) Gff^,M , (2.24) 

where 

Gff,oAu^s) = -e^^ E I in - l,r' \ X,^ \ n,r) ^ (2.25) 

and = — A;TlnX]exp(— /3e„^r) is the grand canonical potential.il The equivalent equations 
for the c-electrons are obtained by just replacing | {n — l,r' \ Xqo- | n,r) p in Eq. ( p.25| ) by 
I (n — 1, r' I Cj^o- I n, r) p. 
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The f-electron GF can be written in the form 

^ m ■ 

Gff,M = - expm Y: , (2.26) 

where uj are the poles and rrij the residues of the GF. There are only eight different uj for 
the f-electron GF, because different transitions have the same energy and the residues of 
some transitions are zero. Each uj = — [Sn-iy — £n,r) corresponds to the lines identified 
with j that appear joining the levels in figure 0, the two lines ui and the single lines M3 and 
u-r represent transitions that are allowed in the absence of hybridization, while the remaining 
ones correspond to transitions that are forbidden in that limit. It is important to notice that 
for a system with given values of E^, Ef and V, the position of the levels in figure § changes 
with the chemical potential /i. In that figure we have Eq = fi and Ef < Eq = fi, and in that 
system the ground state is always the singlet | 2, 9), which has no magnetic moment in the 
absence of field but can have a rather large induced moment because of the proximity of the 
magnetic triplet.!^ 



C. The atomic effective cumulant approximation 

The atomic effective cumulant approximation (AECA) consists in substituting 1^2^/ (z) 
in Eq. (|2.2CI| ) by an approximate M2*^{z) derived from the exact solution of the atomic limit, 
obtained by solving for M2*^(z) in the equation that is the atomic equivalent of Eq. ( p. 20] ). 
One then obtains 

and from the point of view of the cumulant expansion, it contains all the irreducible dia- 
grams that contribute to the exact M2^J {z). It should be emphasized that this diagrams 
contain loops of any size, because there is no excluded site in this expansion, and all the 
filled circles correspond to the same site, although they appear as different vertices in the 
diagram. The difference between the exact and approximate quantities is that different en- 
ergies i?k,(T appear in the c-electron propagators of the effective cumulant M2^J [z), while 
these energies are all equal to Eq in M2*^{z). Although M2^^{z) is for that reason only an 
approximation, it contains all the diagrams that should be present, and one would expect 
that the corresponding GF would have fairly realistic features. 

One still has to decide what value of Eq should be taken. As the most important region 
of the conduction electrons is the Fermi energy, we shall use Eq = n — SEo, leaving the 
freedom of small changes 6E0 to adjust the results to particular situations, but fixing its 
value for a given system when fi has to change to keep the total number of electrons Nt 
fixed, as for example when changing the temperature T. 

Another important point, is that concentrating all the conduction electrons at Eq would 
overestimate their contribution to the effective cumulant, and we shall then reduce the 
hybridization by a coefficient that gives the fraction of c-electrons which contribute most. 
We consider that this is of the order of Vg^, where is the density of states of the free 
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c electrons per site and per spin, and to be more definite we chose irVg^, so the effective 
hybridization constant Va coincides with the usual "mixing strength" A = irV'^g^. This is 
essentially the same choice made by Alascio et.al.0 in their localized description of valence 
fluctuations. Note that Va is only used in the calculation of M2^^{z), and that the full value 
must be substituted in the V that appears explicitly in Eq. ( |2.2CI| ), because the whole band 
of conduction energies is used in G°^{k,z) = —1/ {z — ^ko-)- 

The spectral density p/ = pffla{uj) for several values of T is shown in figure |^ assuming 
the following system's parameters: Ej^^j = Ej = —0.5, /i = 0., a local hybridization V = 0.3 
and a density of states of the unperturbed band electrons given by a rectangular band of 
width 71 centered at the origin. As discussed before, we include adequate universal constants 
into the different parameters, so that all of them are expressed in terms of a single energy 
taken as unit. In the present case we can take that unit to be equal to the total width of 
the band divided into tt. 

The localized energy Ef of the local state is well below the Fermi surface in this figure, 
corresponding to a typical Kondo region. The spectral density pf does not change with T 
below a certain value, which is approximately T = 10^^ in figure ^, because below this T 
the state | n = 2,r = 9) is the only occupied in the atomic limit and therefore M2^^{z) does 
not change by further decreasing T. Other states are occupied at higher T, and the changes 
in M2*^{z) are reflected in p/. The pf obtained by the AECA at the lowest T in the figure is 
basically the same pf of the CHA in the region close to ej = —0.5, but it has also a structure 
close to the Fermi surface (i.e. at = 0) that is absent in the CHA. The spectral density in 
this region has the main characteristics of the Kondo peak, namely its localization and the 
decrease of its intensity when T increases, as can be seen in figure 0. This structure shows 
a pseudo gap, as was obtained using other methods, but with a peak below the Fermi 
energy that is not as sharp and that becomes more complex with increasing T, an effect 
that we believe is due to the use of the atomic model to estimate M|{.'^(z). 

In figure ^ it is plotted the dependence of the spectral density with p showing a main 
structure that follows £/, and a smaller one, corresponding to the Kondo peak, that remains 
fixed at the Fermi surface when the system remains in the Kondo region, i.e. when < 0. 
When Ef approaches the Fermi surface in the intermediate valence region the two parts of 
the spectrum merge, and a single structure then follows ej for > 0. 

The AECA gives a p/ that roughly agrees with the results obtained by other methods,!! 
but the details of the spectral density near the Kondo resonance depend in a very delicate 
way on the behavior of M|{.-^(z) near that region. In the calculation of M2^^{z) all the band 
structure is replaced by a single value Eq, and it is therefore not surprising that the AECA 
results are less precise than those of reference,!! which uses the decoupling of the equations 
of motion followed by the self-consistent determination of the averages resulting from that 
procedure. This last method, as well as the dynamical mean field theory,0 require some 
heavy computation while the AECA is fairly simple from that point of view. As the spectral 
density p/ in the AECA has the same overall behavior shown by the methods mentioned 
above, it can be used to obtain rather reasonable values of many physical properties and of 
their dependence with T and other parameters. The main interest in this method, is that it 
is very natural to make the extension to more complex systems with numerous local states, 
that can be simplified by using the Hubbard operators to project the Hamiltonian to the 
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subspace of states of interest. The use of the AECA would then make it possible to calculate 
properties of those systems without employing too heavy computation. 



III. MAGNETIC SUSCEPTIBILITY AND RESISTIVITY OF FESI 

We discuss here the Kondo insulators within the context of the PAM, employing the 



approximate GF introduced in Section |II Q . We consider the extreme case of U oo, and 
apply it to FeSi, which seems to behave like a typical Kondo insulator.Ha Although this 
compound should be better described with a finite U, we believe that this would not change 
too much the basic properties from the ones presented here, but there would certainly be 
changes in the parameters necessary to reproduce the measured properties of FeSi when a 
finite U is employed. Both the static conductivity cr(T) and magnetic susceptibility xiT) 
vanish when T — 0, and we can describe this dependence of the two properties, as well 
as their initial rapid increase and posterior behavior for increasing T, with the same set of 
values of our model's parameters. The agreement with x{T) is very good at not too high 
T, but to obtain a good agreement at high T one should consider that the parameters are 
affected by thermal expansion in this region. 

To select the total number ritat of electrons per site for our model we first considered 
that in a band calculation of FeSi,Q where all the Fe 3cfAs^ and Si 3s^3p^ states were treated 
as valence electrons, the corresponding 48 valence electrons per unit cell (4FeSi) completely 
filled the lowest 24 valence bands. To make a connection with the much simpler PAM, we 
first notice that a model of two hybridized bands with two electrons per site would also fill 
the lowest band, and that this model was considered to be a good starting point to describe 
FeSi.Ei Therefore it seemed that Utot = 2 would be a good choice for our model. It is easier 
to keep a constant fi to calculate the T dependence of properties, but the corresponding 
changes in ritot where large enough to produce substantial changes in x{T) but minor changes 
in p{T) as the temperature increases. It was then necessary to resort to the more laborious 
calculation at constant ritot- 

The parameters employed in the numerical calculations are numbers that have to be 
scaled later to describe the different particular physical situations, and for brevity we shall 
call them "unsealed parameters" . This unsealed parameters are usually chosen so that one 
physical magnitude, like the band width, is equal to one. In our case we use a rectangular 
band, but fix the unsealed total width 2W = 5 vr, this value (or any other of similar order 
of magnitude) being convenient from the point of view of the numerical calculation. The 
unsealed energy parameters can then be considered to be given in units of 2W/5'n'. In this 
paper we shall present a study of a typical intermediate valence situation with ritot = 2, 
and these two conditions are satisfied by employing values of Ef^^,. and /x close to 4. when 
2W = 5 vr. We shall then use as parameter the common value Ef of Ef ^^ in the absence 
of magnetic field, and fix an unsealed value of Ef = 4. As a final step, the value of fi that 
satisfies Utot = 2 has to be obtained before the x{T) and p{T) are calculated. 

When U —>■ oo, the average of the identity in the space of the local states at site j 

^i;TT + + ^jfio = h , (3-1) 
gives the "completeness" condition 
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(3.2) 



where n^o- and Uf^ are the average occupation number per site of f-electrons states | j,a) 
and I j, 0) respectively. This relation gives the conservation of probability of the local states 
at site j, and it is clear that it should be satisfied before most properties of the system could 
be calculated. The equivalent condition for the uncorrected conduction electrons, obtained 
from the identity Cj^Cj^^ + Cj^^jCj^^ = 1, is satisfied for each a in all the approximations that 
we have employed for the cumulant expansion, but the Eq. ( ^.21 ) is not satisfied for the local 
electrons in the U ^ oo limit .111 The origin of this problem has been already discussed,i3 and 
an exact solution for some families of diagrams, together with a conjecture for the general 
case, has been presented,0 but this treatment is not adequate for the approximate GF we 
employ in this work. We have then applied an "ad hoc" renormalization, multiplying each 
unrenormalized Gjj^(k, z) by a constant x^-, so that the occupation numbers calculated with 
the renormalized Gff^(^(k,z) = x„G^j-^^(k, z) satisfy Eq. ( |3.2| ). The two GF with a = ±1 are 
equal in the absence of magnetic field, so that Eq. ( |3.2| ) is sufficient to determine the two 
identical x^, and in the following section we discuss how to obtain the two different in 
the presence of a magnetic field. 



A. Magnetic susceptibility of FeSi 

The simplest way to calculate the static susceptibility is to consider the one electron GFs 
in the presence of a small magnetic field in the z direction, and obtain rif^cr, nc^a by employing 
Eqs.(|2ri^-^]T^) and the corresponding relations for the conduction electrons. The static 
magnetic susceptibility is then proportional to gf {rif^^ — Uf^-^) + gc (n^j^ — Uc^i) when we 
neglect the Van Vleck contributions, and the calculation of the corresponding approximate 
GFs in the presence of a field is simpler when the gyromagnetic factors gf and gc of the two 
type of electrons are equal. The susceptibility x(T) is then proportional to — divided 
into the magnetic field, where = n + ^c,(t is the total number of electrons per site for 
each spin component a. To calculate the two factors x^ of the "ad hoc" renormalization 
we need of a further relation besides Eq. (|3.2|). We notice that each G//^o-(k, z) would give 
an independent ra/^o? and the two resulting values are in general different in the presence of 
the field, so we ask of them to be equal as the extra condition. The following symmetric 
equations are then obtained 

xt + xin^ = 1 

Di + xi^i = 1 , (3.3) 



where n/^^ and D„ are given in Eqs. (p.l4|j2.19| ) 



To compare the ratio x(T)/x(Tm) with the corresponding experimental value, it is nec- 
essary to chose the energy units to make the fit, and it was found that a convenient method 
was to start the calculations with the given unsealed 2W and Ef and a trial set of the re- 
maining parameters, and then obtain the position T^'-'^ of the calculated maximum of x{T)- 
The constant st defined by Tm = 536 K = st T^^'^ would then give the required scaling, 
so that multiplying all the unsealed energy parameters times sr, would give their absolute 
values in units of degree K. As the maximum of x(T) is rather fiat, it was found that a 
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better procedure was to find an st that would give the best fit of x{T) / xiTm) at moderately 
low T. All the adjustments were made by trial and error and direct comparison of the two 
curves in the plot 

A very good fit was obtained by keeping a constant /i = 4.3 with the 2W = 57r and 
Ef = 4. discussed above, and the relation ritot = 2.0 was approximately satisfied at low T, 
but we could not find a similar set of parameters that would give a good fit for all the values 
of T when requiring a constant Utot = 2.0. As the interval of temperatures is rather large, 
we considered that because of the lattice thermal expansion, some of the parameters could 
also change with T, and for simplicity we employed a linear change in the hybridization 
constant: V{T) = Vq (1 + ay T), where it is sufficient to leave all the scaling in Vq and use 
always the unsealed ay and T in the product ay T. A very good fit shown in figure H was 
then obtained with Vq = 1.8 and ay = —1.2 and a position of the zero-width conduction 
band, employed to calculate M2^^{z), given hy Eq = n — 6Eq with a 6Eq = 0.3 independent 
of T. The scaling of the parameters was given by st = 1750 K ~ 0.1508 eV, so that in 
absolute units 2W = 2.379 eV, Ef = 603.2 meV, Vq = 271.4 meV and 6Eo = 45.24 meV. 

We shall employ the same values of these parameters when comparing the theory with 
the measured resistivity in the following section, so that a consistent description of the two 
properties be obtained. 

It is possible from our calculations to analyze the changes with T of the total occupations 
Uf and ric of the f and c electrons respectively, and we have found that in the calculated 
curve of figure ^ the nj changes from 0.47 at low T to 0.53 at 800 K. These parameter values 
show that the system remains in the intermediate valence region. It is also interesting to 
compare the separate contributions Xf{T) and Xc{T) of the f and c electrons to the total 
xiT) = XfiT) + XciT), and these two quantities are plotted in figure |. Both follow the 
same trend of the total xiT)^ but XciT) is much smaller than Xf(T), so that the choice of 
equal gyromagnetic factors for the two type of electrons should not affect the final result in 
a substantial way. 

B. Resistivity of FeSi 

The dynamic conductivitv o" (u;, T) is related to the current current correlations by the 
well known Kubo formula.^J^l Two-particle GF are then necessary to calculate those corre- 
lations, and to simplify the calculations, Schweitzer and Czychollc3 employed the expression 
of the conductivity for dimension d = cxd as an approximation of the static conductivity for 
d = 3. Only one-particle GFs are then necessary to obtain a (a;, T) in that limit, because the 
vertex corrections cancel out,i3 and we shall follow the same approach. As the hybridization 
is a hopping of electrons between two different bands, it contributes to the current operator 
in the PAM,S but this contribution cancels out in our model because we employ a local 
hybridization Vj,k,cr = V^>- The expression obtained contains explicit sums over k, but it 
is possible to make a further simplification by considering nearest-neighbor hopping in a 
simple cubic lattice,ii@ii and the sums over k can be transformedi^ in integrals over the 
free conduction electron energy e (k). This transformation is possible because in the AECA 
the Gcc,aO<; is a function of k only through the ^(k) = e, as both M|*^(z) and ^-.k.o- = ^> 
are k independent. We can then write 
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pc,a{i^]e) = - lira Im{Gcc,ai'^,uJ + i\'r]\)} , (3.4) 



and the dynamical conductivity is given by 

/"CXD poo 

a {u, T) = Co- / dujf [fri^') - hi^' + t^)] / de pcA^'\ e) pc,a{oji + e) (e) 

UJ J — OD J — OO 

(3.5) 

where Co is a constant discussed below and g1{s) is the density of states of the free conduction 
electrons per site and per spin. 

The static conductivity is then given by 



a{T) = Co r du (- 



dfriuj) 



S{u;) , (3.6) 



duj 
where 

/oo 
ds iPcA^;e)f glie) . (3.7) 
-oo 

The constant 

^ e2 1 2 t2 

Co = TT - - — , 3.8 

n a d 

where a is the lattice parameter of FeSi and t is the hopping parameter of the hypercubic 
lattice that gives the c electron energyll 

^(k) = -(^)Ecos(aM • (3.9) 

We shall generally use a rectangular band with —W < e (k) < W , and to relate t to 
W we consider that the band width y2d2 t of Eq. ( |3.9| ) should be equal to 2W ^ so 
that t = W/\/2d. Replacing d = ?>, 2W = and a = 4.489Ai in Eq. (^) we findil 
Co = 6180/ {n cm). 

If in the AECA of Eq. (|2.21|) we abbreviate Gcc,cr(k, z) = — [a — e(k) + i with real 
a and b (functions of u) defined hj a + i b = z+ \ V {"^ M2^^{z), we find that when 6 — 
then S{uj) ~ O |1/6| if c<j is outside the gap and 5'(ti;) ~ O if u; is inside. This different 
behaviors can give rather different low T limits of a (T) , because the integrand in Eq. ( p. 61 ) 
only contributes in an interval of O |T| around the Fermi energy a; = (our frequency 
variables are given with respect to p). When = is inside the gap we then have a very 
small a (0) if Im M2*g.(0) ~ 0, which corresponds to the case of FeSi at low T, as is discussed 
below. 

To analyze M2*^(z) we notice from Eq. ( p.26| ) that Cj^jrQ^{z) is real when rj = Im [z] — >■ 



except at its only singularities on the real axis, that are the poles at z = uj. It is then clear 
from Eq. ( |2.27| ) that Ml^^^^z) is real on the real axis of the z complex plane,0 except at 
the real solutions oi {u — Eq + p) — \ V p C'fjoA'^) = 0^ where it would have poles with 
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Im M2^^{z) 7^ in their neighborhood. One should then make all the calculations at a finite 
r] and afterwards take 77 — 0. 

In figure ^ we plot the local spectral density of the conduction electrons, namely 

Pc,a{uj) = - lim /m|-^ VG'cc,a(k,u; + i |?7|) I (3.10) 

as well as Im M2^^{uj + i \ri\) as a function of u for the same parameters employed in figure 
but in this plot the variable u, as well as the fixed parameters T = 0.001 and 77 = 0.00001, are 
given in the unsealed energy units. It is clear from figure |^ that for these rj and T the peaks 
of Im M2^^{uj + i \r]\) are very sharp, and that the region where this quantity is appreciably 
large is far away from a; = in units of T. The value of the h defined above is therefore very 
small at that u; = 0, and the conductivity would be extremely low because that value of ui 
is well inside the gap so that S{uj) ~ O If the Fermi surface (u; = ) were inside the 
conduction band, the Siu) ~ O |1/6| and the conductivity would be extremely large, and 
would tend to infinity when ^ 0. The physical reason for this different behaviors at very 
low T is the absence of carriers when = is inside the gap, and the absence of scattering 
mechanisms for the c-electrons when = is inside the band. 

The extreme sharpness of the structure of Im M2^^{z) is a consequence of the atomic 
approximation employed, and to alleviate this character we have added an extra imaginary 
part rja = \ria\ sgn{Im[z]) to its argument : M2*„{z) =^ h/H^^^ioj + irja), so that the poles 
of this quantity become Lorentzians that somehow mimic the effect of the band width. To 
show the effect of this change, the plot of Im M|*^(z) is also shown in figure |^ for the same 
parameters used above but for an unsealed rj + rja = 0.001. The Im M2^^{z) has now an 
appreciable value at = and reaches inside the conduction band, and cr(T) is much larger 
than before, but it is still rather small because T is small and only an interval of order of T 
contributes to cr(T). 

Addition of rja to the argument of M2^^{z) leads to similar effects as those already obtained 
by Mutou and HirashimaBa through "introducing a small imaginary part T to the conduction 
electrons", i.e. replacing z = iujhjz + iT sgnioj) in the GFs G//^o-(k, 2;) and Gcc,ai}^iz). 
Their justification is the existence in real systems of scattering processes due to phonons and 
impurities, and we should also consider this mechanisms as contributing to the irja. Within 
this interpretation one could also consider a temperature dependence of r/a, but we have not 
implemented this change in the present calculation. 

It seems clear that the basic scattering mechanism in our calculation of the PAM's 
a (T) is the hybridization, because the otherwise free conduction electrons are scattered by 
the localized f electrons through this interaction. This is apparent if we notice that the 
relaxation effects are described by the imaginary part of the usual self-energy Scc,o-(k,z), 
defined through 

Gec,a(k,2) = -{^-£(k)-S,e,.(k,2)}-^ , (3.11) 

and that the exact relation Ecc,a(k,;z) = - | V^(k) p M|{^(k,;z) follows from Eq. The 
relaxation mechanism of the c-electrons is then provided by the hybridization, and the self 
energy is independent of k in the AECA: Tjcc,u{z) = — \ V M2*^(z). 
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As discussed above, the vanishing resistivity obtained in our approximation when n is 
inside the conduction band and T — is caused by the atomic character of our effective 
cumulant M2^^, and the introduction of a finite rja moderates this effect. In our calculation 
of p{T) = l/a{T) we have chosen rja to obtain a reasonable agreement with the measured 
values for T above 50 K, because the conduction by ionized impurities would dominate 
<y(T) at lower temperatures.!! In figure |^ we plot p(T) for the same parameters employed 
to calculate xiT) in figure ^ , but for an unsealed rja = 0.0008. The plot is compared to the 
experimental values measured! by Schlesinger et.al., which we obtained by digitalization of 
their plot. We have preferred these values to those of other authorsi^'0 because the higher 
values of p(T) seem to indicate a better sample with less impurities. The agreement is fairly 
good above 50 K, and the much lower experimental resistivity below this T is attributed to 
the ionic conduction, as mentioned above. The impurity scattering is T independent and 
could be considered to be included in the value of rja we used, but the phonon scattering 
would increase with T, and one should start with a smaller rja and add a contribution that 
has the T dependence corresponding to this process. This effect could explain the difference 
between the two curves above T = 200 K. 

We point out that the typical values of rja we have used would not affect the values of 
x{T), although they were essential in fitting p(T) to the experimental values. 



IV. CONCLUSIONS 

In this work we have used approximateil one-electron GFs of the PAM in the limit U ^ oo 
to study the static magnetic susceptibility x{T) and electrical resistivity p(T) of FeSi, which 
has the typical behavior of the Kondo insulators. The total number of electrons per site 
ritot = 2 we have selected seems a good choice for FeSi, and at each T it is then necessary 
to find the value of the chemical potential p that gives this ritot- 

The exact GFs of the PAM are formally derived from their cumulant expansion,^! and are 
expressed in terms of an effective cumulant M2^J (k,z) , that is given by the contribution of all 
the proper diagrams of the exact GF Gfj^^ik, z). To obtain approximate GFs, we substitute 
M2^/ (k, z), whose calculation is as difficult as that of G ff^^^ik, z), by the corresponding 
quantity of an exactly soluble model, which is the same PAM but with a conduction band 
of zero width. The use of this atomic limit makes the corresponding effective cumulant 
^2*cr(^) independent of k, and the approximate GFs have the same structure obtained by 
the dynamical mean field theory.E^I While this last method uses self-consistency conditions 
to obtain the final GFs, we use physical considerations to chose the parameters that define 
M2^^{z). The approximation we employ does not give results as accurate as that theory, but 
has the advantage that the numerical calculation is fairly rapid, and that it can be extended 
to more complex systems with a large number of local states, by projecting the large space 
of these states into a manageable subspace of interest with the use of Hubbard operators. 

We have shown that the spectral density p/ of the approximate GF employed in this 
work has, at low T and in the Kondo region, a structure on the Fermi surface that can 
be interpreted as the Kondo peak, a feature that was missing in our cumulant expansions 
when only cumulants up to fourth order were employed. This behavior is essential for the 
calculation of the system's properties, and as the spectral densities obtained by the AECA 
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have the overall behavior that is expected from the model,El the use of this approximation 
would then be expected to give, without too heavy computation, rather reasonable values 
of many physical properties and of their dependence with T and other parameters. 

These approximate one-electron GFs of the PAM in the limit U ^ oo, are used here 
to describe the static magnetic susceptibility x(T) and electrical resistivity p{T) of FeSi, a 
compound which has the typical behavior of the Kondo insulators. We have selected a total 
number of electrons per site ritot = 2, because it seems a good choice for FeSi, and at each 
T it was then necessary to find the value of the chemical potential /i that gives this ritot- 
We tried to describe the system employing a typical intermediate valence situation, and this 
was achieved with a rectangular band of total width 2W = Svr and an unperturbed energy 
Ef = 4. of the f-electrons, requiring a chemical potential /i ~ 4 to satisfy ritot = 2. These 
unsealed values correspond to measuring the energy parameters in units of 2W/5tt. 

The susceptibility x(T) is proportional to the difference —n^ between the occupation 
number of up and down electrons, and we have fitted the calculated x(^)/x(^m) to the 
experimental results,0 where x{Tm) is the maximum value of the susceptibility. We have 
found that the contribution of the conduction electrons to x{T) is much smaller than that 
of the local electrons for the employed parameters (cf. figure and the assumption made 
of equal gyromagnetic factors for the two type of electrons has therefore no major relevance. 

To find the factor st that transforms the unsealed parameters into absolute units, we 
have scaled the calculated x(T)/x(Tm) curve so that it agrees with the experimental one at 
low T. It was not possible to find a set of parameters that would give a good agreement 
at both low and high T, and we considered that because of thermal expansion some of the 
parameters could change with T. Assuming a linear dependence V{T) = Vq (1 + ay T) of 
the hybridization parameter V it was then possible to obtain a rather good agreement in the 
whole T range (cf. figure |^). Although this is a very crude model of these effects and one 
should also expect changes of Ef, it shows that the dependence of the model's parameters 
with thermal expansion should not be neglected. 

To calculate p(T) we used an expression valid for infinite dimensionj^l and assumed a 
nearest neighbor hopping in a simple cubic lattice to transform the sums over k into integrals 
over the unperturbed energies £:(k) of the conduction electrons. As a further non-essential 
simplification we used a rectangular band for the corresponding spectral density g^. 

We have shown that the relaxation mechanism of the c-electrons is provided by the 
hybridization, and the most likely situation in the AECA is that the conductivity vanishes 
as T — > when /i is inside the gap and tends to infinity when /i is inside the band. Here 
we show that this result is a consequence of the atomic character of the AECA, that makes 
M2^^{z) imaginary only in the neighborhood of its finite number of poles that are usually 
far apart from the Fermi surface {uj = 0). This behavior is moderated by introducing 
an extra imaginary part rja to its argument (cf. figure J3), whose effect is similar to that of 
scattering of the c-electrons by phonons and impurities.a The resistivity at low temperatures 
is probably determined by ionic conduction, and if we choose rja so that the calculated 
resistivity coincides with the experimental one close to T ~ 50K, a very good agreement is 
obtained up to T ~ 2007^. The departures at higher T could be attributed to the phonon 
scattering. 

In summary, we have shown that the spectral densities obtained by the AECA have 
the overall behavior that is expected from the PAM, and we can therefore expect that 
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this approximation would give, without too heavy computation, rather reasonable values of 
many physical properties and of their dependence with T and other parameters. We have 
employed the AECA to study FeSi assuming a typical intermediate valence situation of the 
PAM in the C/ — > oo hmit, and we were able to give a good description of both the static 
resistivity and magnetic susceptibility of FeSi as a function of T employing a single set of 
the model's parameters. 
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FIGURES 



FIG. 1. Typical cumulant diagrams for one-particle GF, where the filled and empty circles 
represent the f-electron and c-electron free propagators respectively, and the lines joining them 
represent the hybridization interaction, (a) The diagrams of the chain approximation (CHA) for 
the f-electrons, represented by the filled square to the right, (b) As (a) but for the c-electrons, 
represented by an empty square, (c) A more complicated diagram, with cumulants of fourth and 
sixth order. 

FIG. 2. The energies e^.r = -E'n.r — of the twelve eigenstates | n, r) of the atomic limit are 
represented in this figure, and those corresponding to different occupations n = 0, 1, 2, 3 are drawn 
in different columns. The index r that characterize the states is written above the correspond- 
ing levels, and the lines joining different levels are identified by numbers i, showing the possible 
transitions Uj that contribute to the GF. As with Eqs. ( |2.^j2.9| ), the frequencies have the chemical 
potential /i subtracted, so that the Fermi surface corresponds to w = 0. 

FIG. 3. The spectral density of the f-electrons obtained with the AECA for several values 

of T, employing z = lo + if] with tj = 10~^. The system has the following parameters: Ef = —0.5, 
/i = 0., r = 0.001, a local hybridization V = 0.3 and a density of states of the unperturbed 
band electrons given by a rectangular density of states of width vr centered at the origin, all in 
the same energy units. The effective cumulant was calculated employing a reduced hybridization 
Va = A = V'^ and a frequency with a small imaginary part r]a = 0.003 (further details of the 
introduction of rja are given in Subsection 1IIB| ). 



FIG. 4. The spectral density of the f-electrons for the same system of figure |3| but with fixed 
T = 0.001 and for several values of the chemical potential /U. It shows the crossover from the Kondo 
region to the intermediate valence region. 

FIG. 5. The reduced susceptibility x{T)/x{Tm) as a function of T, where x{Tm) is the maximum 
value of x(^) • The open circles correspond to the experimental values, taken from reference |^. 
The full line corresponds to the model's calculation for a rectangular band with a total width 
2W = Stt, an f electron energy of Ef = 4., a total number of electrons ntot = 2.0, a T dependent 
hybridization V{T) = 1.8(1. — 1.2T) and the zero-width conduction band employed in the model 
located at Eq = /j, — 6Eq with a T-independent 5Eo = 0.3. These unsealed values have to be 
multiplied by the scale factor st = 1750 K = 150.8 meV to express them in absolute units: 
2W ~ 2.379 eV, Ef = 603.2 meV, Vq = 271.4 meV and 6Eo = 45.24 meV. 

FIG. 6. The separate susceptibilities Xf{T) (full line) and Xc{T) (dotted line) of the f and c 
electrons respectively, as a function of T and for the same parameters of figure ^. 
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FIG. 7. The spectral density of the conduction electrons Pcai^^) times vr (full line) and the 
imaginary part of the approximate effective cumulant M2^^(z) (dotted and dash-dotted lines)are 
plotted in a logarithmic scale as a function of u; for T = 0.001, r] = 0.00001, rja = 0. (full and 
dotted line), 1] + r]a = 0.001 (dash-dotted line) and the same parameters of figure |5[ The values 
of LJ, T, r] and r]a are given in unsealed energy units, and have to be multiplied by the scale factor 
St = 1750 K= 150.8 meV to express them in absolute units. 

FIG. 8. The resistivity ,(T) in Q cm, measured by Schlesinger ed.ahi (triangles) and calculated 
(full line) with Eq. ( |3.6D and the same parameters employed in figure ^, but with rja = 0.0008. 
The values of uj and r] are given in unsealed units, and have to be multiplied by the scale factor 
St = 1750 K= 150.8 meV to express them in absolute units. 
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TABLES 
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TABLE I. The properties of the twelve eigenstates | ri, r) of H. are given. The columns are 
labeled by the number of electrons n, the name r of the state, the z spin component Sz of the state 
and the value of en,r = -£'n,r ~ ^/^j where is the energy of the state | n,r). 
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